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SECTION  1 
INTRODUCTION 


BACKGROUND 

The  Air  Force  is  interested  in  determining  the  impact  of  airport  activities  on 
local  ambient  air  quality.  Currently,  an  air-quality  computer  model  such  as 
the  Air  Quality  Assessment  Model  (AQAM)  is  used  to  evaluate  the  air-quality 
impact  of  airports,  both  in  their  current  configurations  and  with  proposed  op- 
erational or  other  changes.  In  most  air-quality  models,  including  AQAM,  a 
Gaussian  dispersion  formulation  which  is  particularly  adaptable  to  small  dis- 
tances and  pollutant  travel  times  associated  with  airports  is  used.  Models 
based  on  the  superimposi tion  and  addition  of  Gaussian  plumes  from  point,  ’Ine, 
or  area  sources  have  rather  difficult  output,  which  must  be  interpreted  for 
assessment  of  environmental  impact.  The  direct  computer  output  of  AQAM,  for 
instance,  are  tables  with  up  to  312  X-  and  Y-coordinates  of  pollution  concen- 
trations for  each  of  five  pollutants.  These  tables  are  both  time  consuming 
to  read  and  meaningless  to  someone  unfamiliar  with  air-pollution  analysis. 

A contour-plotting  routine  is  therefore  required  to  graphically  produce  lines 
of  equal  concentrations.  These  contour  plots  can  be  used  in  conjunction  with 
a map  of  an  Air  Force  Base  and  vicinity  for  air-pollution  dispersion  analysis. 


OBJECTIVE 

To  minimize  the  efforts  of  AQAM  users  in  obtaining  computer-contoured  displays 
of  AQAM  predictions,  the  contouring  pa;kage  should  be  fully  interfaced  with 
AQAM.  This  could  be  accomplished  by  a contouring  input  tape  generated  during 
runs  of  AQAM.  Thus,  the  purpose  of  this  research  effort  was  to  develop  and 


r 


ip 


implement  a plot  input  tape  mechanism  to  interface  the  AQAM  short-  and  long- 
term models  (refs.  1 and  2)  with  the  contouring  plot  package  developed  at  the 
Civil  Engineering  Research  Facility  (CERF)  of  the  University  of  New  Mexico  and 
to  test  the  contour  plot  package. 

SCOPE 

A software  package  was  designed  to  accept  data  from  cards  or  from  a contour 
plot  input  tape  generated  on  an  AQAM  run.  The  user  should  have  control 
through  an  input  command  structure,  over  which  data  will  be  plotted,  and  the 
specific  contour  levels,  their  tensions  and  dashline  patterns  that  would  be 
used  in  each  contour  plot.  It  was  decided  that  a driver  program  which  would 
interface  the  user's  intentions  --  expressed  in  terms  of  input  instructions, 
user  data,  and  the  contouring  code  --  and  the  system  plot  package  was  prefer- 
able to  an  approach  in  which  the  user  would  have  to  code  a driver  program  for 
each  different  application  and  recompile  the  contouring  program. 

The  data-management  mechanism  for  AQAM  calculations  is  developed  in  section  2 
of  this  report.  Since  the  input  data  matrix  of  predictions  is  typically 
equal  to  17  by  17,  it  was  necessary  to  provide  linear  bivariate  interpolation 
of  the  original  AQAM  output  to  effectively  double  or  triple  the  dimensions  of 
the  input  matrix.  This  capability,  expansion,  is  also  presented.  The  con- 
touring algorithm  is  described  in  section  3.  Two  other  important  code  fea- 
tures for  contour-curve  fitting  with  splines  under  tension  and  a dashline- 
patterned  contour-line  drawing  capability  are  also  discussed.  Graphic  display 
capabilities  and  several  examples  are  provided  in  section  4;  detailed  informa- 
tion on  the  required  input  instructions  for  the  AQAM  contour-plotting  package 
is  given  in  section  5;  and  conclusions  and  recommendations  are  contained  in 
section  6.  Examples  of  input  instructions  are  also  provided  (appendix  A) 
followed  by  examples  of  AQAM  calculations  for  Scott  Air  Force  Base  (appendix  B). 


1.  Wangen,  L.E.,  and  Rote,  D.M.,  A Generalized  Air  Quality  Assessment.  Model 
for  Air  Force  Op, 'rations,  AFWL-TR-74-54 , Air  Force  Weapons  Laboratory, 
Kirtland  Air  Force  Base,  New  Mexico,  May  1975. 

2.  Wangen,  L.E.,  and  Rote,  D.M.,  A Generalized  Air  Quality  Assessment.  Model 
for  Air  Force  Operations  - An  Op, rater's  Guide,  AFWL-TR-74-304 , Air  Force 
Weapons  Laboratory,  Kirtland  Air  Force  Base,  New  Mexico,  February  1975. 
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SECTION  2 

DATA  MANAGEMENT  OF  AQAM  PREDICTIONS 


It  was  necessary  to  develop  an  AQAM  plot  restart  tape  format  which  would  be 
well  suited  to  the  generalized  model  (ref.  3).  AQAM-predicted  concentration 
data  fall  into  two  categories,  grid  concentration  data  and  special  receptor 
concentration  data,  as  shown  below: 


Grid 

Concentration 

Data 

Special  Receptor 
Concentration  Data 

Pu 

P.2 

P.3  

••  P 1 n 

CSR, 

P2  1 

P22 

P23  

..  P2n 

CSR? 

P31 

P32 

P 3 3 

..  P3n 

CSR, 

pk, 

Pk2 

pkl 

••  Pkn 

CSRe 

The  mapping  of  grid  concentration  data  into  a list  of  concentrations  and  the 
concentration  of  special  receptor  data  are  as  follows: 

C.  is  the  concentration  of  pollutant  for  receptor  grid  cells  j = 1,  m and  m = 
” t h 

17  x 17;  CSR.J  is  the  concentration  of  the  i special  receptor  at  and 

l = 15. 

(Cl,  C2,  C3  ....  Cm,  CSR i , CSR?  ...  CSR.) 

The  relationship  between  and  (shown  above)  is  as  follows: 

Ci  C2  ...  C 1 7 C 1 8 ...  Cm 


Pii  P i 2 . . . P i n P?i  ...  P^h 

where  m=17xl7,N=17byK=17.  Thus,  P^  is  mapped  into  C^  where  n = 1, 
N and  k = 1,  K and  j = k + (n  - 1)*K. 

3.  Menicucci,  David  F . , Air  Quality  Assessment  W<  iel  (AQAM)  ita  faction 

and  Operator's  Guile  (Graft),  AFWL-TR-75-307 , Air  Force  Weapons  Labora- 
tory, Kirtland  Air  Force  Base,  New  Mexico,  December  1975. 
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(Time  Block  V (Time  Block? 
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|T  ime  Block  N 


(a)  Restart  Tape 


Header 

Pol lutant 

Pollutant 

Pol  lutant 

Pol lutant 

Pol lutant 

Record 

HC  Data 

_ J 

NO  Data 

Particulate 

X 

Data 

(h) 

Time  Block 

Parameter 

Forma  t 

Variable 

T i me 

(no) 

I CT IME 

Month 

(AID) 

I MONTH 

Weekday  - 0,  Weekend  1 

(no) 

1WD0F 

Title 

(A10) 

ITITLE(G) 

Number  of  Pollutants 

(no) 

HPOL 

Pol lutants 

X-  and  Y-Coordinates  of 
Lower  Left  Corner  of 

(A10) 

I POL  L ( 5 ) 

Receptor  Grid 

(F10. 3) 

XYCRG(2) 

Number  of  Rows 

(no) 

NROW 

Number  of  Columns 

(no) 

NCOL 

X-  and  Y-Grid  Spacing 

(F10. 3) 

XYSP 

Humber  of  Special  Receptors 
X-  and  Y-Coordinates  of 

(no) 

NSRP 

Special  Receptor 
Meteorological  Data 

(Flo. 3) 

XYSPE ( 30) 

Stabi 1 i ty  Category 

(no) 

ISC 

Wind  Speed  (m/sec) 

(F10.3) 

ws 

Wind  Direction  (deg) 

(Flo. 3) 

WD 

Temperature  (rT) 

(F10.3) 

TEMP 

Mixing  Depth  (m) 

(F10.3) 

DMID 

(c)  Header  Record 


1 Concentration 

1 Concentration  1 

Concentration 

Concentration 

Data  for 

Data  for 

Data  for 

Data  for  Total 

Aircraft 

Airbase 

Environment 

(Combined) 

(d)  Pollutant  Data  Block 


Figure  1.  Structure  of  AQAM  Contouring  Input  Tape 
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The  contouring  input  tape  contains  several  time  blocks  as  shown  in  figure  la. 
Each  time  block  contains  a header  record  and  a pollutant  data  block  for  each 
of  five  pollutant  types  as  illustrated  in  figure  lb.  The  format  of  the  header 
record  is  shown  in  figure  lc.  Each  pollutant  data  block  contains  concentra- 
tion data  for  the  four  source  classifications:  aircraft,  airbase,  environ- 

ment, and  total  (fig.  Id). 

The  AQAM  receptor  grid  represents  a sampling  domain  in  two  dimensions,  in 
which  model  calculations  are  made.  The  lower  left-hand  corner  of  the  grid, 
the  number  of  rows  and  columns,  and  the  spacing  between  mesh  points  are  chosen 
by  the  user.  The  typical  AQAM  grid  size  is  17  by  17,  but  rectangular  grids 
can  be  used  as  long  as  the  total  receptors  are  less  than  312.  The  smoothness 
of  the  contours  is  dependent  both  on  the  tension  of  the  contour  and  the  rank 
of  dimension  of  the  AQAM  receptor  grid.  The  expansion  capability  allows  the 
user  to  double  or  triple  the  dimension  of  the  original  receptor  grid  through 
a simple  linear  bivariate  interpolation  algorithm.  It  is  generally  felt  that 
the  interpolated  matrix  contains  no  more  real  information  than  did  the  orig- 
inal matrix.  AQAM  predictions  at  special  receptor  locations  within  the  recep- 
tor grid  are  likely  to  show  much  higher  concentrations  than  predictions  at 
mesh  points  because  of  their  proximity  to  sources  when  taking  into  account 
meteorological  conditions.  The  expansion  capability  is  beneficial  when  one 
must  use  a receptor  grid  that  is  small  (less  than  12  by  12)  or  when  there  are 
spikes  in  the  grid  predictions. 
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SECTION  3 

CONTOURING  ALGORITHM 


The  development  of  the  AQAM  contouring  program  involved  the  establishment  of  a 
data-preparation  mechanism  and  computer  software,  modification  and  implemen- 
tation of  existing  software  packages,  and  concatenation  of  all  software  with 
system  plot  packages.  A literature  search  was  performed  but  very  little 
useful  information  on  computer  contouring  was  found.  The  contouring  programs 
available  from  computer  software  companies  were  either  too  expensive,  too 
complicated  in  general,  produced  contour  plots  of  insufficient  quality  and 
accuracy,  or  required  too  much  effort  in  terms  of  input  on  the  part  of  the 
applications  programmers. 

CONREC  CONTOURING  CODE 

Communications  with  Mr.  Tom  Wright  of  the  National  Center  for  Atmospheric 
Research  (NCAR)  indicated  that  a contouring  code  named  CONREC,  in  use  at  NCAR, 
might  be  adapted  for  use  in  conjunction  with  AQAM.  CONREC,  authored  by  Mr. 

Tom  Wright,  is  a Fortran  callable  subroutine  for  two-dimensional  contour 
plotting  (ref.  4).  CONREC  accepts  a two-dimensional  array  of  input  data  to  be 
contoured.  For  each  contour  level  the  algorithm  searches  along  the  perimeters 
of  the  input  matrix,  interrogating  successive  mesh  points.  In  the  event  the 
mesh  point  to  the  left  is  less  than  or  equal  to  the  contour  value  which  is 
also  greater  than  the  mesh  point  to  the  right,  the  algorithm  finds  the  star- 
ting position  for  a given  contour  level.  A subprogram  is  called  to  trace  the 
contour  from  this  starting  position  by  linearly  interpolating  the  neighboring 
mesh  points.  The  inequality  between  left  and  nght  mesh  points  and  book- 
keeping in  the  subprogram  which  traces  the  contour  line  guarantee  that  no  * 

contour  will  be  drawn  twice.  The  internal  rows  of  mesh  points  of  the  input 
array  are  similarly  searched  for  starting  positions,  and  each  contour  is 
traced  in  turn  from  its  starting  position.  This  procedure  is  repeated  for 
each  contour  level.  The  contour  lines  which  are  to  be  labeled  are  drawn  first. 


4.  Wright,  T.J.,  Utility  Plotting  Pmjnms  at  NCAR , Atmospheric  Technology, 
National  Center  for  Atmospheric  Research,  September  1973. 
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In  the  smoothing  version  of  this  algorithm,  the  linearly  interpolated  points 
for  each  contour  line  are  smoothed  by  splines  under  tension  (ref.  5).  The 
smoothed  contour  lines  are  then  labeled  by  a software  dashline  subroutine. 


The  operating  system  and  the  Fortran  compiler  at  NCAR,  however,  are  not 
standard  CDC  software.  Computer  codes  at  NCAR  are  not  compatible  with  the 
standard  CDC  6600  compilers  like  those  at  AFWL.  More  specifically,  nonstan- 
dard input,  output,  and  shifting  operations  are  handled  differently  by  the 
Fortran  compilers  at  these  installations.  It  was  therefore  necessary  to 
modify  the  NCAR  routines  which  had  NCAR  system  dependencies  and  add  codes 
appropriate  to  the  standard  CDC  Fortran  FTN  4.2  compiler.  Linking  the  con- 
touring code  to  the  AFWL  system  plot  package  was  accomplished  with  a set  of 
stop-gap  routines. 

CONREC  was  also  modified  to  allow  for  greater  user  control,  and  hence  flexi- 
bility, in  the  type  of  contours  generated.  Data  structures  which  allow  the 
user  to  specify  unequal  contour  levels,  a tension  parameter  used  in  smoothing, 
and  the  dashline  pattern  for  each  contour  line  were  implemented. 

Smooth  contours  which  do  not  overlap  must  be  obtained  by  both  adjusting  the 
tension  parameter  and  by  increasing  the  resolution  of  the  input  array.  A 
subroutine  which  increases  the  input  mesh  data  from  N by  K to  2N  - 1 by  2K  - 1 
was  programmed  by  averaging  neighboring  points  for  bivariate  linear  inter- 
polation. Other  resolution-increasing  techniques  such  as  bicubic  spline 
interpolation  tend  to  produce  oscillatory  interpolation  for  typical  AQAM 
concentration  data  because  of  the  rapidly  changing  derivatives  in  the  concen- 
tration surface  and  therefore  were  not  considered. 

SPLINES  UNDER  TENSION 

Two  important  requirements  for  contour  maps  are  that  contour  lines  not  inter- 
sect and  each  contour  line  be  smooth.  A number  of  curve-fitting  techniques 
whi<~h  allow  for  curve  smoothing  have  been  developed,  but  when  they  are  applied 


5.  Cline,  A.K.,  "Scalar-and-Planar-Val ued  Curve  Fitting  Using  Splines 
Under  Tension",  Comm  ACM,  Vol.  17,  No. 4,  April  1974,  pp.  218-223. 
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to  contour  plotting,  smooth  contour  lines  which  overlap  or  intersect  result. 
A.K.  Cline  (ref.  5)  using  splines  under  tension  has  developed  techniques  for 
planar  curve  fitting.  Both  opened  (nonperiodic)  and  closed  (periodic)  curves 
are  considered.  The  important  characteristic  of  splines  under  tension,  as 
opposed  to  other  curve-smoothing  and  interpolation  techniques,  is  that  with 
sufficient  tension  a curve  can  be  drawn  as  close  as  desired  to  polygonal  lines 
drawn  through  the  knots  and  still  be  smooth.  A unique  feature  of  the  AQAM 
contour-plotting  program  is  the  capability  of  specifying  the  tension  parameter 
for  fitting  smooth  contour  lines  for  a particular  contour  level.  This  is 
important  in  contouring  air-pollution  concentrations.  Typically,  there  are 
large  spikes  in  the  concentration  surface.  When  contours  are  drawn  at  levels 
which  are  changing  rapidly--corresponding  to  large  partial  derivatives  in  the 
concentration  surface--the  resolution  of  tie  data  for  that  contour  is  typi- 
cally very  low.  This  means  that  the  contour  curve  is  based  on  only  a few 
knots.  An  increase  in  resolution  may  provide  more  knots  on  which  to  base  the 
contour.  To  circumvent  the  need  for  unreasonably  high  resolution  data, 
splines  under  tension  can  be  used  effectively  to  insure  that  contours  which 
are  based  on  a small  number  of  knots  do  not  intersect.  In  terms  of  AQAM 
contouring  applications,  the  implication  is  that  low  tension  can  be  specified 
for  low  contour  levels  and  greater  tension  can  be  specified  for  higher  contour 
levels  where  there  is  a possibility  that  the  low  resolution  of  the  data  to- 
gether with  insufficient  tension  will  produce  an  intersection  of  contour 
lines.  The  dimensionless  tension  parameter  ranges  from  1 to  30.  Figure  2 
illustrates  the  effect  of  tensions  1 to  9 on  curves  drawn  through  four  knots. 


DASHLINE  PATTERNS 

Some  choices  of  contour  levels  inevitably  produce  small  contours  without 
labels  which  indicate  the  height  of  the  contour.  These  may  be  level  curves 
near  the  top  of  a spike  or  local  maxima.  Dashline  patterns  can  be  used 
effectively  when  this  situation  occurs;  i.e.,  the  pattern  of  the  line  identi- 
fies the  level  even  though  it  has  no  label.  Thus,  the  user  by  specifying  a 
dashline  pattern  for  each  contour  level  can  enhance  the  clarity  of  the  contour 
plot.  These  patterns  and  their  numbers  are  given  in  figure  3. 
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fENSN-6.  0 


1ENSN-5. 0 


TENSh  4.0 


TENSN-9.0 


1EN5N-7. 0 


TEN5N-8.  0 


Figure  2.  Effect  of  Tension  Parameter 


Figure  3.  Dashline  Patterns 
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SECTION  4 

GRAPHIC  DISPLAY  TECHNIQUES 


In  the  interest  of  obtaining  clarity  in  the  presentation  of  AQAM  calculations 
an  investigation  was  made  into  other  computer  graphic  display  techniques. 
Figure  4 illustrates  the  superimposition  of  line-drawn  contour  plots  on  half 
tone  plots  of  the  same  data.  A halftone  graphic  used  in  conjunction  with 
line-drawn  contours  enhances  the  visual  recognition  of  the  data  being  dis- 
played. 

Halftone  shading  can  be  especially  useful  when  the  data  are  very  noisy.  AQAM 
predictions  on  a receptor  grid  are  generally  somewhat  noisy;  spikes  in  the 
concentration  surface  can  occur  because  of  the  magnitude  of  a source  and  its 
proximity  to  a grid  point.  These  spikes  create  some  difficulty  in  choosing 
contour  levels.  Without  a careful  choice  of  unequal  contour  levels,  these 
spikes  make  reading  a contour  plot  rather  difficult.  By  using  the  halftone 
shading  technique,  the  spikes  can  easily  be  seen  even  if  they  are  very  narrow 

Another  graphic  display  of  two-dimensional  arrays  of  data  is  the  line-drawn 
surface  plot  (fig.  5).  Hidden-line  elimination  allows  the  surface  to  be  dis- 
played as  though  it  were  opaque  or  solid.  Figure  5 also  demonstrates  the 
effect  of  increasing  the  resolution  (from  17  x 17  to  33  x 33)  by  the  linear, 
resolution-increasing,  expansion  technique. 
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Figure  4.  Halftone  and  Contour  Lines 


(b)  33  x 33  Resolution 


Figure  5.  Line-Drawn  Surface  Plot 
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SECTION  5 
INPUT  INSTRUCTIONS 


The  input  cards  are  used  as  follows:  to  retrieve  a specified  time  block  from 

the  restart  data  tape;  to  select  options  such  as  various  sources,  pollutants, 
and  special  receptors;  to  perform  special  operations  such  as  different  contour 
levels  and  tensions;  and  to  alter  the  content  of  the  title  as  given  on  the 
input  tape.  Command  cards  are  utilized  to  allow  as  much  flexibility  as  pos- 
sible. They  are  described  below: 


GETTIME  DEFAULT  - search  for  a specified  time  block  and  process  with 
default  options. 

GETTIME  CARD  - read  data  on  punched  cards  immediately  following  the 
command  card.  This  is  used  when  the  restart  data  tape 
is  not  available. 

TITLE  - read  text  on  the  next  punched  card  to  override  the 

title  given  by  the  restart  tape. 

OPTION  - read  the  following  punched  cards  which  contain  the 

data  on  contour  levels,  tensions,  and  dashline  pat- 
terns. This  option  is  effective  only  during  the  cur- 
rent GETTIME. 


END  OF  OPTION 
END  OF  JOB 
NO  LEGEND 


- end  of  contour  option  data. 

- end  job. 

- eliminate  the  legend  from  the  plot. 


The  minimum  input  data  cards  for  a successful  run  are  any  number  of  GETTIME 
cards  and  an  END  OF  JOB  card.  Between  10  and  30  contour  levels  are  generated, 
and  the  number  of  levels  is  dependent  upon  the  input  data. 

GETTIME  DEFAULT  Card  (one  card  required) 


Columns 

Variable 

Format 

Description 

1-7 

GETTIME 

(A10) 

Input  Command 

11-17 

DEFAULT 

(A10) 

Input  Command 

21-22 

LS 

(A2) 

ST  - Short-Term  Model 
LT  = Long-Term  Model 

23-26 

ICTME 

(A4) 

Time  Block  Descriptor 
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Columns 

31-34 


41-45 


51-55 


Variable 

ISO 


IPO 


IEXP 


Format 

(411) 


(511) 


(511) 


Description 

Source 


Column 

Source 

IS0(1)  31 

Aircraft 

I SO ( 2 ) 32 

Airport 

I SO ( 3 ) 33 

Environment 

ISO (4)  34 

Iota  1 

If  I SO ( n ) =1  , the  n 

sources  wi 1 1 

be  plotted  and  all 

others  will 

be  ignored. 

Default:  all  sources  will  be 

plotted . 

Pollutant 

Column 

Source 

IP0(1)  41 

CO 

I PO ( 2 ) 42 

HC 

I PO ( 3 ) 43 

NO 

x 

IPO (4)  44 

Particulate 

IPO (5)  45 

S0? 

If  IP0(n)  = l , the  n 

^ pollutants 

will  be  plotted  and  all  others 

will  be  ignored. 

Default:  all  pollutants  will  be 

plotted . 

Expansion  of  Matrix 

Column 

Source 

I EXP ( 1 ) 51 

CO 

I EXP ( 2 ) 52 

HC 

I EXP ( 3 ) 53 

NO 

X 

I EXP ( 4 ) 54 

Particulate 

I EXP ( 5 ) 55 

SO 

2 

If  I EXP ( n ) =0 , the 

original  matrix 

with  rank  NCOL,  NROW  will  be 

contoured . 

If  I EXP ( n ) = 1 , the 

matrix  wi 1 1 be 

expanded  to  (2  x NCOL ) - 1 , 

(2  x NROW) -1 . 

If  IEXP(n)=2,  the 

matrix  from 

the  first  expansion  will  be 

expanded  again. 

i .e.  , NCOL = 1 0 

NR0W=1 2 

lEXP(n)=0 

gives  10  x 12 

IEXP(n)=l 

gives  19  x 23 

IEXP(n)=2 

gives  37  x 45 

No  defaul t values . 
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Columns  Variable  Format  Description 


Columns  Variable  Format  Description 

61-65  I SPEC  (511)  Special  Receptor 


Column  Pollutant 

ISPEC(l)  61  CO 

I SPEC ( 2 ) 62  HC 

I SPEC ( 3 ) 63  N0x 

r SPEC ( 4 ) 64  Particulate 

1SPEC  5)  65  SO 

? 

If  I SPEC ( n ) = 1 , special  receptors 

of  the  n1'*1  pollutants  will  be 
plotted . 

Default:  Special  receptors  will 

be  ignored. 


GETT1ME  CARD 

Card  (minimum  of 

three  cards  required) 

Card  1 

Columns 

Variable 

Format  Description 

1-7 

GETTIME 

Command 

11-14 

CARD 

Command 

Card  2 and  successive  cards  necessary  for  this  group  of  data  are  to  be  read 
in  as  namelist  group.  This  means  that  keywords  are  required  to  be  punched 
for  corresponding  values.  See  reference  manual  (ref.  6).  The  following  key- 
words are  used  in  this  group  and  are  explained  in  the  restart  tape  format: 


ICT1ME 

- 

not  necessary 

IM0NTH 

- 

octal 

IWD0E 

- 

integer 

ITITLE(8) 

- 

octal 

IP0LL ( 5) 

- 

octal 

XYCRG( 2) 

- 

real 

NR0W 

- 

integer 

NC0L 

- 

integer 

XYSP 

- 

rea  1 

NSRP 

- 

integer 

ISC 

- 

integer 

WS 

- 

real 

WD 

- 

real 

TEMP 

- 

real 

DMID 

- 

real 

I SORCS ( 4 ) 

- 

octal 

NPOL 

- 

integer 

6.  Fortran  Fxtrnde!  V.-rr.f.  n ! !u  Control  Data  Corporation, 

Cyber  70  Computer  Systems  Models  72,  73,  74,  76,  7600  Computer  System 
and  6000  Computer  System,  1974. 


ISO 

- integer 

I POL 

- integer 

IEXP 

- integer 

ISPEC 

- integer 

Cards  numbered  from  n + 1 on  are  the  concentration  data. 


TITLE  Card  (two  cards  required) 


Columns 

Variabl e 

Format 

Description 

1-5 

TITLE 

Command 

1-80 

INTIT 

(8A10) 

Command 

CONTOUR  OPTION 

Card 

Card  1 

Columns 

Variable 

Format 

Description 

1-6 

OPTION 

Command 

Card  2 

Columns 

Variable 

Format 

Description 

11-20 

CONL 

(F10.3) 

Contour  Level 

21-31 

TENS 

(F10.3) 

Tensions 

31-41 

LNP 

(HO) 

Dashl ine  Pattern 

Card  3 is  the  same  s CONTOUR  OPTION  Card  2 for  next  contour  level. 


Card  n-1  (Note:  Contour  levels  must  be  read  in  ascending  order.  A maximum 

of  20  may  be  inputted.) 

Card  n 

Columns  Variable  Format  Description 

1-10  END  OF  OPT  Command 


NO  LEGEND  Card  (one  card  op‘ iona 1 J 
Columns  Variable  Format 

1-9  NO  LEGEND 


Description 

Command 


IB 


END  OF  JOB  Card  (one  card  required) 

Columns  Variable  Format  Description 

1-10  END  OF  JOB  Command 

As  described  under  GETTIME  CARD,  some  of  the  variables  being  read  in  as 
a namelist  group  must  be  converted  from  characters  to  display  code  as  follows 


A 

01 

P 

16 

4 

B 

02 

Q 

17 

5 

32 

C 

03 

R 

18 

6 

33 

D 

04 

S 

19 

7 

34 

E 

05 

T 

20 

8 

35 

F 

06 

U 

21 

9 

36 

G 

07 

V 

22 

+ 

37 

H 

08 

W 

23 

- 

38 

I 

09 

X 

24 

★ 

39 

J 

10 

Y 

25 

/ 

40 

K 

11 

Z 

26 

( 

41 

L 

12 

0 

27 

) 

42 

M 

13 

1 

28 

Blank 

45 

N 

14 

2 

29 

Comma 

46 

0 

15 

3 

30 

Period 

47 

Consider  as  an  example  that  one  wishes  to  have  the  month  on  the  plot  read 
SEPT  for  September.  The  namelist  input  will  then  read: 

1M0NTH=1 9051 6204 04 5454 54545B . 
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SECTION  6 

CONCLUSIONS  AND  RECOMMENDATIONS 


(1)  The  software  package  is  extremely  fast;  it  requires  only  about  3 
sec  of  execution  time  per  contour  plot  for  a 17  x 17  matrix. 

(2)  The  input  command  structure  requires  a minimum  amount  of  input  informa- 
tion but  allows  the  user  to  directly  control  the  display  of  AQAM  cal- 
culations. 

(3)  The  contouring  package  is  fully  interfaced  with  the  AQAM  via  the  con- 
touring input  tape  generated  during  AQAM  runs. 

(4)  Plots  produced  by  the  contouring  package  can  be  automatically  labeled, 
if  desired,  with  appropriate  information. 

(5)  The  package  has  been  implemented  on  the  CDC  6600  Computer  System  at 
Eglin  Air  Force  Base  and  is  connected  via  intercom  to  a DATA  100  sys- 
tem at  Tyndall  Air  Force  Base.  Plotting  will  be  done  on  a Calcomp  936 
plotter  at  the  Air  Force  Civil  Engineering  Center,  Tyndall  Air  Force  Base. 

(6)  The  AQAM  contouring  package  represents  the  state-of-the-art  in  computer 
graphic  display  of  two-dimensional  contours. 

(7)  The  user  can  easily  control  the  content  of  the  graphic  output. 

Recommendations  for  future  research  and  development  which  are  beyond  the 

scope  of  this  technical  effort  are  as  follows: 

(1)  The  halftone  technique  described  in  this  report  should  be  imple- 
mented . 

(2)  An  effort  should  be  made  to  develop  methods  for  analysis  of  graphic 
display  of  concentration  data.  An  example  of  this  type  of  analysis 
is  the  computation  of  areas  enclosed  by  contours  for  each  contour 
level.  Another  might  be  the  computation  and  display  of  the  rate  of 
change  of  the  concentrations  over  the  grid. 

(3)  Techniques  for  the  overlay  of  digital  cartography  of  an  airbase  map 
and  vicinity  on  contour  plots  should  be  explored. 

(4)  Techniques  for  producing  camera-ready,  report-quality  graphics  at 
minimum  expense  should  be  developed. 

(5)  The  contouring  code  should  be  adjusted  so  that  very  small  contours 
are  not  plotted. 
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APPENDIX  A 

EXAMPLES  OF  INPUT  INSTRUCTIONS  FOR  CONTOURING  TAPE  AND  CARD  INPUT 


r 


CONTOURING  WITH  TAPE  OR  DISK  FILE  INPUT 

The  sample  contour  plots  were  generated  from  predictions  made  by  AQAM  for 
Scott  Air  Force  Base,  Illinois.  Figure  A1  demonstrates  the  result  of  GETTIME 
DEFAULT  with  automatic  contouring;  figure  A2  demonstrates  the  use  of  GETTIME 
DEFAULT  with  contour  OPTIONS  for  the  same  data.  The  Scott  graphics  in  appendix 


B illustrate  the  technique  for  the  overlay  of  a line-drawn  contour  plot  on  a 


halftone  graphic  of  an  airbase  map. 

The  input  instructions  for  contouring 

with 

tape  input  (figs.  A1 

and  A2)  are  as  fol lows; 

1 

7/8/9 

2 

G 0 

ST0600 

1 1 

3 

TITLE 

4 

SCOTT  AFB  AND  AFCS 

AND  C130 

(BASE  AND  VICINITY)  GRID  SPACING  .6  KM 

5 

G 0 

ST0600 

1 1 

6 

TITLE 

7 

SCOTT  AFB  AND  AFCS 

AND  Cl 30 

(BASE  AND  VICINITY)  GRID  SPACING  .6  KM 

& 

OPTION 

9 

100. 

2.5 

0 

10 

250. 

2.5 

4 

11 

500. 

5.5 

4 

12 

900. 

5.5 

4 

13 

END  OF  OPT 

14 

END  OF  JOB 

15 

6/  7 / 8/ 9 

The 

first  GETTIME  DEFAULT 

results  in 

automatic  contour  level  selection  by  the 

program.  Contour  options  for  the  same  data  are  specified  in  the  second  GETTIME 
DEFAULT. 


Dashline  patterns  can  be  used  to  readily  identify  contour  levels  of  special 
interest.  In  the  examples,  the  higher  contour  levels  are  drawn  with  solid 
lines.  The  relationship  of  the  airbase  map  and  vicinity  to  the  contour  plot 
is  determined  by  the  initial  choice  of  the  location  of  the  AQAM  (X,Y)  grid. 
Photoreproduction  techniques  ensure  an  accurate  overlay.  The  airbase  map  is 
piiotoreduced  and  screened  before  being  overlaid  with  the  contour  plot. 
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CONTOURING  WITH  CARD  INPUT 


The  following  example  illustrates  the  full  input  instruction  set  for  card 


input. 

The  contour  plot  produced  by  these  input  instructions  would 

be  the 

same  as 

that  shown  in  figure  A2. 

1 

7/5/9 

45 

4 . 5 7 8 £ <•  C 2 

2 

**  O 

6 

2.147E*Ql 

3 

{INCAR 

7 

5.  liitE  + u 1 

4 

ICTIME=2T,  3T  , 27,  ?7,33,c7,34,27,?7,45. 

<.  *1 

J.h97^;.  1 

5 

INONTH=lr>f"c  (Jvi  1 3»C2»C5tl8  i‘t3i 

7 9 

1 • -•  6 ir  «•  j ^ 

6 

5l 

1 . r 2 

7 

*VCR'T=244.?4J  , <.266.5';  1 , 

51 

2.<-37E*C2 

1 

N°OW= 1 7 , 

7? 

2 . 4 5 8 c.  *■  2 

i 

NCGL-17, 

7 3 

2 .ub7cK  2 

L 

*YSD=. 6u , 

C 4 

,^4CE<-C  2 

1 1 

NSPP=9 , 

55 

2.137c*? 2 

12 

”30=6, 

5 o 

1 3 

WS=1. C l , 

57 

1 .5l7E-'6 

1 <♦ 

W0  = 1T-  . -o , 

5 3 

3.21 flE-„4 

1 5 

T£RP=62. 1 , 

59 

1 * -*r  2E  -C  2 

16 

OH  13=  1.  j."1) 

6 * 

2 . 732E-C  1 

1 7 

TE  xp( i»  , 

5 1 

0.49  7E*Cr 

1 1 

TSP€~(0»=0* 

c o 

1 .1  9 8E+  r 1 

19 

IS0(*> =1 , 

63 

5.667E*01 

TO0(l > = 1 , 

5 *■ 

A « .,  * 5c  ♦ & 1 

j « 
U ▲ 

NPOLl-5 1 

6 5 

1 . 3 9 C E ♦ 0 2 

<.2 

o 5 

1 . 4 7 4c  * rl  2 

2 3 

1 .4316-26 

t>7 

2.420E+T2 

24 

1 .527F-C4 

4 a 

2.5 5lE*o2 

25 

5 ,26  0E-P  7 

6 9 

2.tn68E*J  2 

2o 

9.570E-C2 

7 u 

2 . 5 1 9 E ♦ u 2 

? 7 

742E-U 

71 

2.c67E*G2 

2 * 

3.7-»2E»JC 

72 

1 .96?E»P2 

29 

1.276C»C 1 

7 3 

3" 

■».2l7c*0l 

**  *♦ 

i.-37E-u6 

71 

6.399E*?1 

75 

j . 23  4E-C  4 

32 

l«I54t*,w2 

7 

2.693r-p2 

3 3 

1 .495E* j 2 

7 7 

^ .52  BE-C  1 

3 4 

1 . 2 

7 3 

- . °8  1 E 5 r 

35 

7.:56E»n2 

79 

2 ,_5CE*.  1 

'6 

2.2ai£*T2 

a r 

b.317E+ii 

?.29<.E  + ,,2 

5 1 

1 .■’69E*'-? 

3 * 

?.193E*32 

*2 

1.984E+C2 

*9 

1 3 

?.564£f u2 

•»  0 

1 • jk  3»” C 6 

6 t 

8 76E4-C  z 

4 1 

2 • 1 4 6E  - 0 <♦ 

3 5 

Ml  5E*  C 2 

4t 

a .467c-C 3 

* 6 

’ ./?CE*u2 

*.3 

1 • 3 0 w c. " L 1 

a 7 

? .-C9EO 2 

4U 

1 .721E»lp 

3 3 

2 • l **5£  ♦ „ 2 

89 

1 . 6 (5  2 E ♦ n 2 

24 


r 


9C 

i : -> 

16* 

1 . " 39E**  7 

91 

1. 2? it-” 6 

1 ’ 5 

2.224£ *c  2 

13  1 

6,1o7E*C2 

92 

3 • 4 1 ftc  “ j 4 

137 

1 . 5 5 4 £ * ] 2 

13? 

-.926E*(,2 

1’ 

5.52 7E-* 2 

1 _ 3 

1.  * 9 4-1*97 

1°  3 

3 . 1 1 2r * l 2 

9-+ 

1 . 1 9 2 E * u . 

1 ' 9 

3. 5 j 9c*  0 1 

1*4 

1 .925r*C£ 

95 

1 .r.  5 4E*]  1 

*4  . 

2 1 r * * 1 

1*0 

1.237c*C2 

3 6 

..726c*r  1 

1-1 

1 ) 6 

7 . 9 4 6 E ♦ : 1 

9 7 

1 .r  A6c*!;£ 

1 ^ 2 

1 . '1 6£-C  7 

1«  7 

5.w54r*0i 

j 3 

1 .9 44c* 3 2 

1.  3 

1 . 4 3 c;  _ - r 2 

16  6 

T . 1 7 3E*  0 1 

99 

2 . 7 1 C 5 ♦ l 2 

1.4 

1 . '» 9 « t * G 5 

139 

1 .985E*Q 1 

1L  0 

7 . 1 4 2 c * 0 2 

1-5 

7 . hi  5c*  . 1 

19t 

1 .?4lE*n  1 

1.  1 

7 . 1 9 1 E ♦ C 2 

1 -6 

1. 7u3E*C  2 

1 91 

7.79  6c*li: 

1'  2 

7. 953c*  j 2 

1 .7 

2 . ^ 6-:* : 2 

19  c 

1 ' o 

?.5b2E*J2 

1 U 3 

-».2ob£*r2 

1 u 3 

1 -» 

° .121E *C  c 

1 . 9 

4 .5156*02 

194 

2 . 37QE-0  4 

1 « 5 

1 .o95E*72 

lr  ? 

7.  U 36*i  . 

1 9 c 

j.7HE*0l 

1.  o 

i . 32  : c.*r  2 

151 

3 . . 3 1 c ♦ : 2 

1 9 6 

1 • L 5 *£  * ' t 

1*  7 

152 

2.21 6E * 2 2 

197 

3.9716*52 

1„  3 

lotfil-.’ 

15  3 

1 .9o7£*C2 

1 > 3 

o.613c*C2 

1L  9 

i .^24c-: 3 

1 5 4 

1 . L 3 4E  ♦:  7 

15? 

?.7o6E*"2 

XI  f* 

1 .272£-C  1 

155 

7.4i0^*:  i 

?"  r 

1 . 6 3 3 £ * C Z 

ill 

2.343E*».^ 

156 

> . .53c*.  1 

ci  2 

1 .w41c*12 

112 

2.766E* jl 

15  7 

3 . ‘.43c*C  1 

2.  c 

o.77SE*C 1 

il  3 

•-UlM 

15  <3 

?C  3 

* . 1 4 6£  * J 1 

114 

1 .(J2<>E*',  2 

159 

1 . >j  2E-C  5 

2.  4 
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Figure  A1 . Automatic  Contour 
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Figure  A2.  Contour  with  Options 
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